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NOMENCLATURE 


All of the data is referenced to fuselage body axes according 
to right-handed sign conventions. 


Si 

b 

CG 

Cl 

c. 

c 

c 

d 


f 

g 

h 

h 


Horizontal distance from pivot axis to spring line of action, 
feet 

Normal acceleration, g 

Longitudinal acceleration, g 

Reference span, feet 

Center of gravity, fraction of chord 

Lift coefficient 

Pitching moment coefficient 

Reference chord, feet 

Viscous retarding moment 

Perpendicular distance from the CG to the oscillation axis, 
feet 

Load applied to tail tiedov«m, Ibf 
Frequency of oscillation, cycles per second 
Acceleration of gravity, ft/sec^ 

Altitude, feet 

Vertical, component of the perpendicular distance from the CG 
to the oscillation axis, feet 

Moment of inertia about the axis of rotation, slug-ft^ 

Moment of inertia about roll axis, slug-ft^ 

Moment of inertia about pitch axis, slug-ft^ 
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Ij Moment of inertia about yaw axis, slug-ft^ 

k Spring constant, Ibf/ft 

Ka Flow amplification factor for angle of attack, 1.1 

MAC Mean aerodynamic chord 

MSL Mean sea level 

p Roll rate, de/sec 

q Pitch rate, deg/sec 

^ Dynamic pressure, Ibf/ft^ 

r Yaw rate, deg/ sec 

5 Reference area, ft^ 

Tj. Thrust coefficient 

u Control vector 

V Velocity, ft/sec 

W Weight, Ibf 

Wg Aircraft and supporting system weight, Ibf 

X State vector 

X Amplitude of oscillation 

Horizontal distance from to feet 

Distance of alpha vain and an accelerometer forward of CG, 

Xa feet 

n 

z Measured observation vector 

1 Vertical CG position from support pivot axis, feet 

Z Vertical distance from the support pivot axis to the aircraft 

reference line, feet 

a Angle of attack, deg 

6 Angle of sideslip, deg 

5 Logarithmic decrement 


V 1 


Sq Elevator deflection, deg 

9 Pitch angle, deg 

5 Damping ratio 

Roll attitude, deg 

Natural and damped frequencies of oscillation, rad/sec 

Subscripts 
m Measured 

n Number of cycles 

q Rotary derivative, per rad 

a Static derivative, per deg 

6g Control derivative 

0 Bias 

Tg Test equipment 
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INTRODUCTION 


Reliable estimates of stability and control derivatives are 
needed in the design of new aircraft to expand the flight envelope by 
flight testing and for research of control systems and aircraft 
handling qualities. Confidence in the values of the derivatives is 
enhanced when there is agreement between values obtained from flight 
tests and wind txmnel tests. 

In order to obtain reliable and accurate values of the sta- 
bility and control derivatives, the Dryden Flight Research Center 
CDFRC) of NASA has developed a technique for extracting the derivatives 
from flight data. This technique is implemented by a set of Fortran 
computer programs (reference 1) that is based on a modified maximum 
likelihood estimator that uses the Newton- Raphson algorithm to perform 
the required minimization of the derivatives. 

DFRC has an ongoing flight test program in which the deriv- 
atives are extracted from flight data of selected aircraft . As a part 
of the program, a PA-30, a light twin-engine general aviation aircraft, 
has been flight tested and the stability and control derivatives 
obtained using the DFRC method. The derivatives were compared with 
those obtained from wind tunnel tests of PA- 30 aircraft in the NASA 
Langley full scale tunnel. The comparison revealed significant dif- 
ferences of derivatives determined by the two methods (see figures 5 
through T), The differences were primarily in the values of the 
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longitudinal stability derivatives. Since specific and detailed docu- 
mentation of the instrumentation used in the flight tests was not 
recorded, it was difficult, if not impossible, to cite a reason for 
the differences in the values of the derivatives. To resolve the dis- 
crepancies it was necessary to repeat the flight tests. 

The second flight tests were conducted as a joint venture 
between the DFRC and California Polytechnic State University as a part 
of the Graduate Student Program. For these tests, particular attention 
was given to the calibration and documentation of the instrumentation 
used. Also, a spring oscillation method (reference 2 ) was used to 
determine more accurate values of the moments of inertia. The same 
aircraft was used for both the Langley wind tunnel (reference 3) and 
flight tests. The test conditions were primarily the same. When con- 
ditions varied, corrections were made in order to make the comparisons 

valid . 

The purpose of this report is to present the results of the 
new flight test program. Data was obtained with the aircraft in zero, 
half, and full flap configurations in level unaccelerated flight with 
the landing gear retracted. The flight tests also included flight with 
the conditions of zero flaps and zero thrust. The data was analyzed 
using the modified maximum likelihood technique to extract the longi- 
tudinal stability and control derivatives. The derivatives were 
plotted as functions of angle of attack using various graphical 
arrangements to show variations of wind tunnel and flight determined 
values at zero flap settings. Also, data was displayed to show the 
effects of flap deflection and thrust variation on the longitudinal 
stability derivatives. 
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DESCRIPTION OF AIRCRAFT AND INSTRUMENTATION 

The aircraft tested was a light twin-engine low-wing monoplane 
with retractable landing gear. A three view drawing of the PA-30 is 
shown in figure 1, and the reference axes is shown in figure 2. The 
geometric and mass parameters of the aircraft are listed in Tables 1 
and 2. The wing airfoil section is a NACA 642 A215- airfoil . The wing 
has five degrees of dihedral with zero twist and is placed at two 
degrees of incidence with respect to the fuselage reference line. The 
horizontal tail is an all -moveable type, with an overall deflection 
range from 16 degrees trailing edge up to 5 degrees trailing edge down. 
The horizontal tail has a trailing edge trim tab that moves in the same 
direction as the tail. The tab has a deflection ratio (tab deflection 
to tail deflection) of approximately 1.5. 

The only modification to the aircraft was the addition of a 
hydraulic control system and a remote pilot/control system used in 
other flight test efforts. 

The instrumentation consisted of a standard package used to 
measure stability and control parameters . This package included 
three-axis angular rate gyros, attitude gyros, and linear acceler- 
ometers in addition to various air data instruments. The locations 
of the instruments requiring positional correction before the data is 
used in the computer analysis are listed in Table 3. 

The auigle of attack and angle of sideslip data were measured 
using vanes mounted on an instrumentation boom located approximately 
3.3 and 3.64 feet, respectively, forward of the right wing tip. 
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The dynamic pressure (q) , velocity, and altitude were calcu- 
lated, using a Fortran computer program, from the static and dynamic 
pressures sampled by the pitot static system onboard the aircraft. 

The pitot tube probe was located approximately 14 feet from the 
aircraft centerline under the left wing. The two static ports were 
located on the sides of the rear fuselage. A radar altimeter, with a 
range of from zero to 5000 feet, was also used to check the calculated 
altitude. 

The control deflections of the stabilator, ailerons, and rudder 
were also recorded along with left throttle position, engine RPM, and 
manifold pressures. The instrumentation analog information was routed 
through signal conditioning circuits and a pulse code modulator (PCM) , 
and subsequently transmitted to ground based recording equipment by way 
of an FM/FM transmitter. All of the instruments were calibrated approx- 
imately one week prior to the flight tests. 

FLIGHT MANEUVERS 

All the data was collected on two flights that yielded a total 
of 97 longitudinal flight maneuvers. Three maneuvers were neglected 
because of excessive time dropouts and data spikes in the output. 

Three more maneuvers were dropped for exceeding tolerable angle of 
attack bounds. The remaining 91 maneuvers were used for the analysis. 

The longitudinal stability and control maneuvers consisted of 
a set of five standard elevator pulses of approximately seven seconds 
duration. The magnitude of the elevator impulse was such that it 
yielded a total change of angle of attack from 3 to 6 degrees. Tnese 
maneuvers were followed by two series of a ”3-2-l-l”-multistep elevator 
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pulse of approximately 11 seconds duration. Both the standard elevator 
pulse and the ”3-2-l-l"-multistep elevator pulse are illustrated in 
figure 3, which includes an example of an in-flight elevator input (5^) 
of each type. 

The "3-2-l-l"-raultistep elevator pulse technique is a maneuver 
that fulfills certain frequency requiranents needed in the parameter 
identification technique used. Theoretically this impulse contains 
frequencies both above and below the aircraft’s natural frequency of 
the short-period mode (see reference 4). 

The maneuvers were performed with the four flight conditions or 
aircraft configurations most often used within the aircraft's flight 
envelope. These configurations included settings of zero flaps, half 
flap (15“), and full flaps (27") at a nominative thrust coefficient, 
and a setting of zero flaps with a thrust coefficient of zero. A set 
of maneuvers was also performed at approximately 300 feet over a flat 
level surface (a dry lake bed) to compare the computed pressure 
altitude with the altitude measured by the radar altimeter. 

The maneuvers were flown in non-turbulent conditions at 
altitudes from 7500 to 13000 feet MSL. Data was obtained at several 
angles of attack for each configuration with the aircraft flying 
straight and level. The angle of attack ranged from .5" to 8“ for the 
zero flap with nominative thrust coefficient flight condition. For the 
one half and full flap settings with nominative thrust coefficient 
flight condition, angles of attack ranged from 3.5" to 8.5". The 
angles of attack for the flight condition of zero flap setting and zero 
thrust coefficient were 1.0", 5.5", and 10.0". 
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All of the flight maneuvers were flown with the landing gear 
retracted, and without the use of stability augmentation equipment. 

The flight tests were conducted with a center of gravity position of 
16% (±2%) MAC, but was corrected to a reference position of 10% MAC 
for analysis . The weight ranged from 3600 to 3450 pounds during the 
two flights, and was known to be within 2% to 3% for each maneuver. 

DATA PROCESSING 

A number of different methods and techniques were used to 
manage and maintain the quality of the flight data, and every attempt 
was made to maintain the integrity of the outputs. 

The recorded flight parameters, which are listed in Table 4, 
were uniformly passed through a 40 hertz passive filter pn_bpard the 
aircraft. The data was sampled by a 10-bit pulse code modulation 
(PCM) digital system at 200 samples per second. The sampling rate was 
reduced to 100 samples per second and an analytic computer program was 
used to filter the buffet noise from the data channels requiring it. 
The filter used a notch at a frequency of 19.4 hertz and a third order 
low-pass at 20 hertz. All of the known phase shifts and spikes in the 
data that were caused by instrumentation or sampling anomalies were 
corrected or removed. The data was ultimately thinned to 25 samples 
per second for analysis. 

The air data which included the parameters of true velocity, 
geopotential pressure altitude, corrected dynamic pressure, and 
corrected static pressure were calculated from the measured static and 
dynamic pressures . 
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The moments of inertia (I^, ly. I^) were calculated using the 
data obtained from the spring oscillation tests of the aircraft. The 
experimental method is shown in Appendix B. The numerical results are 
in Table 2. 


METHOD OF ANALYSIS 

An iterative method of coefficient variance minimization was 
used to find the longitudinal stability and control derivatives from 
flight data. This iterative technique, sometimes called the Newton- 
Raphson algorithm, was incorporated into a general parameter 
estimation computer program. 

The mathematical formulation of this method can be described 
in probabilistic terms. For each possible estimate of the unknown 
parameters, a probability can be obtained that the calculated aircraft 
response time histories will take bn those values actually observed. 
Those estimates are chosen, and the equation's coefficients are 
adjusted so that the probability is maximized. The describing 
equation's coefficients are the stability and control derivatives. 

The best maximization comes when the initial probable values are 
closest to the actual observed values . This process is referred to as 
a maximum likelihood formulation of the problem. 

The maximum likelihood estimation technique can be made to 
include a priori information that can come from wind tunnel studies, 
previous flight tests, and other sources of predicted derivatives. The 
method that uses a priori information is called the modified maximum 
likelihood estimator (MMLE) and is explained in detail in reference 1. 
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In this analysis of the flight test data a priori information was not 
used. 

The full-scale wind tunnel obtained stability derivatives were 
used as starting values for the computer analysis. A complete 
description of the digital computer program (Ntt’lLES) used in the 
extraction of the derivatives is in reference 5. A brief outline of 
the equations of motion used in f^LE3 is included in Appendix A. 

The maximum likelihood estimator also contains an objective 
measure of the validity of the estimates of the stability and control 
derivatives. These uncertainty levels are proportional to the Cramer- 
Rao bounds described in reference 6. These bounds are analogous to the 
standard deviations of the estimated stability derivatives. The 
greater the uncertainty level of the coefficient, the greater the 
uncertainty will be of the estimated derivative. The best uncertainty 
approximation comes by comparing confidence levels of the same deriv- 
atives obtained from different flight maneuvers. Thus by comparing 
uncertainty levels of derivatives obtained from several maneuvers, the 
validity of the stability and control derivatives can be checked. 

A set of stability and control derivatives as a function of 
angle of attack are chosen by fairing the derivative values from 
plotted results at a particular angle. The fairings teike into account 
the derivative grouping and the uncertainty levels which lead to a 
choice of a single representative derivative value at the corresponding 
angle of attack. The faired derivative values from the wind tunnel and 
previous flight tests are listed in Tables 5 and 6, respectively. 
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RESULTS AND DISCUSSION 


The nondiraensional stability derivatives 
and determined from the flight data, and the full-scale Langley 

wind tunnel results were plotted as functions of angle of attack 
(alpha). These results along with representative, faired, values from 
the previous flight tests are graphically superimposed for easy com- 
parison. The graphs clearly illustrate an improvement in the results. 

The uncertainty levels of each derivative are indicated on the 
graphs as a vertical bar. The magnitude of the vertical bars is mul- 
tiplied by five to improve viewing . 

A typical example of a pulse and "3-2-1-1" longitudinal 
maneuver time history of angle of attack, pitch rate, pitch attitude, 
normal acceleration and elevator input is shown in figure 4. The 
dotted lines in these plots are the maximum likelihood method curve 
fit to the actual time history. As can be seen, there is little dif- 
ference between the two curves . 


Zero Flap Deflection, 

Constant Thrust Setting 

The thrust settings for these maneuvers were chosen by the 
pilot to maintain the desired angle of attack. 

The negative value of Cjj, shown in figure 5a indicates a stable 
aircraft. The graph also shows that the magnitude of increases at 
larger values of angle of attack for the flight test data, but 
decreases for the wind tunnel data. The results of the second flight 
test are shown to.be nearer to the values of the wind tunnel results 
than are the results of the first flight test. The C(^ vs alpha curves 
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shovm in figure 5b indicates good agreement between the derivatives 
obtained from the second flight and wind tunnel tests. The 
determined from the first flight test differed as much as 27% lower 
than the tunnel derivatives. The results of the second flight tests 
are clearly an improvement. 

The longitudinal rotary derivative, in figure 5c remained 

negative, constant, and stable over its angle of attack range. There 
were no wind tunnel estimates available for comparison for . How- 

q 

ever, the second flight test results are shown to be consistently 

larger in magnitude than the first test results. 

The vs alpha curve shown in figure 5d reveals that both 

°e 

the first and second flights had curves with negative slopes. The 

values of the derivatives from the second flight are shown to be nearer 

to the values of the wind ttmnel estimates. As shown in figure 5e, the 

values of Ct very closely agreed with the wind tunnel estimates for 

the range of angle of attack tested. Good agreement between the values 

of Ct . determined from wind tunnel and the first flight was obtained 
oe 

only near zero alpha. Excellent agreement for values of for the 

®e 

entire range of alpha tested is shown for wind tunnel and second flight 
determined values. In fact, the most significant improvement of second 
flight results over first flight results is obtained in the values of 



The trimmed elevator position, corrected to the 

reference CG position of 10% MAC, is plotted as a function of angle of 
attack in figure 5f. The pitching moment coefficient, C„,^, and lift 
coefficient, C^ , biases are plotted as a function of angle of attack 
in figures 5g and 5h, respectively. These biases are calculated values 
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in the pitching moment and lift coefficient equations in the maximum 
likelihood procedure. 


Flap Deflection Effects 

Figures 6a through 6h show the effects of flap deflection on 
the stability derivatives. The derivative values from the second 
flight test for half and full flaps are graphically superimposed with 
the zero flap derivative curve. The half and full flap derivative 
curves from the first flight were also graphically superimposed for 
easy comparison. There were no wind tunnel values available for the 
half and full flap conditions. 

The second flight showed more variation between the half and 
full flap derivative values of and than did the first flight. 

As can be seen in the and curves, figures 6a and 6b, the half 
flap derivative values did not change significantly from the no flap 
values. On the other hand, the for the full flap condition were 
greater than those for the full range of alpha, and magnitudes of 
for the full flap condition were lower. 

Both the Cij,^ and derivatives become more negative in mag- 

nitude at lower angles of attack as flap deflection increases from zero 
to half, but does not change as the deflection increases from half to 
full position. See figures 6c and 6d. 

Figure 6e reveals that flap position does not significantly 

affect Cl values at low alpha. However, there are appreciable dif- 
5e 

ferences in the values of Cf „ determined from first and second flight 

*'6e 


tests. 


Figure 6f shows the effect of flap deflection on which 

®trim fl3p deflection increases. The values of 


shows a decrease of 6 
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and Cl^ are also plotted, figures 6g and 6h, to show flap 
deflection effects. 

* 

Thrust Effects 

The effects on the longitudinal stability and control deriv- 
atives caused by varying thrust are depicted in figures 7a through 7h. 
Each graph shows the derivatives obtained for the second flight at T^. 
of approximately zero, along with the first flight and wind tunnel 
curves at zero T<,. The nominative thrust coefficient from the second 
flight is thought to be the same as the wind tunnel value of 0.1. Both 
the second flight nominative and wind tunnel thrust coefficient curves 
®re graphically superimposed to show the change in derivatives. 

Figure 7a shows that the values of 0^^ obtained from both flight 
tests are less negative than the wind tunnel values for zero T^.. There 
is no significant difference in the values for nominative and zero 
T(. conditions. Therefore, the conclusion is that thrust variation does 
not significantly affect the value of for the range of T^ tested. 
Figure 7b reveals that for zero T^, obtained from the second flight 
test are closer to the zero wind tunnel values for small angles of 
attack than are values from the first flight test. The opposite is the 
case for higher angles of attack. For the second flight test the 
nominative T^. derivatives are larger than the zero T,, values. This is 
also true for the wind tunnel nominative and zero T(. values. The con- 
clusion supported by this data is that increase of thrust causes an 
increase of Cj^. 

As seen in figure 7c, values determined from the second 
flight test are less negative for the condition of zero T^ than for 
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nominative T^.. Values of Cn determined from the first flight test are 

even less negative. Again, there are no wind tunnel estimates for 

As shown in figure 7d, for both the wind tunnel and second 

flight tests, values of Cjjjj. for the nominative condition are lower 

Oe 

than for the zero condition. First flight tests values of for 

e 

the zero Tq condition are larger than either the wind tunnel or second 
flight test values. The conclusion drawn from this data is that an 
increase of thrust decreases the values of . 

There is very little difference in the derivative shown 

in figure 7e, as T^ changes. This is also the case for the wind tunnel 
estimates. The conclusion is that thrust variation has little effect 


on C 


L<Se- 


The trimmed elevator deflection is shown for the zero con- 
dition in figure 7f. The moment coefficient and lift coefficient 
biases are shown in figures 7g and 7h. Both biases did not change 
appreciably as the thrust coefficient changed. 


CONCLUDING REMARKS 


A complete set of longitudinal .stability and control deriv- 
atives were obtained . from 91 maneuvers. The conditions of flight that 
were analyzed included zero, half, and full flap configurations with a 
nominal thrust coefficient, and with no flaps and zero thrust coef- 
ficient, The flight itime histories of all the maneuvers had good 
matches with the modified maximum likelihood estimator (MMLE3) pre- 
dictions, which subsequently led to good derivative estimates. Knowing 
the aircraft’s vertical and longitudinal center of gravity, to within 
4$ and 2% respectively, was a major factor in obtaining good matching 
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of the time histories. The second flight estimates agreed with wind 
tunnel estimates much better than those of the first flight. The e.xact 
cause for the differences between the first and second flight results 
is not known. 

Specifically, the conclusions drawn from this study are: 

1. The.Cjj^, C|^, C , and estimates at the zero flap con- 

*^e °e 

dition had closer agreement with the wind tunnel estimates than did the 
first flight estimates. The most significant improvement of the second 
flight results over the first was in . 

2. An increase of flaps from zero to half position shows little 

effect on and , but causes Cjjj and to decrease at lower 

angles of attack. 

3. An increase of flaps above half position causes to increase 
and Cj^ to decrease. 

4. Thrust variation has little effect on C_ , Ci . , C^i , and Ci 

'“a *^Oe o *^0 

for the range of T^, tested. 

5. An increase of thrust causes an increase in Ci , but a decrease 

in Cm and C„ . 

"q '"oe 



TABLE 1 


GEOMETRIC CHARACTERISTICS OF THE PA- 30 AIRPLANE 


WING: 

Area, including aileron and flaps, ft^ 178 

Span, ft 35.98 

Mean aerodynamic chord, ft 5 

Dihedral, leading edge, deg. 5 

Incidence, deg. 2 

Length of flap, each, ft 9.2 

Length of aileron, each, ft 6.3 

Total aileron area, ft* 14.1 

Total flap area, ft* 20.2 

STABILATOR: 

Total area, ft* 52.5 

Span, overall, ft 12.5 

Mean aerodynamic chord, ft 2.7 

VERTICAL TAIL; 

Fin area, ft* 9 

Rudder area, ft* 5.2 

Rudder mean aerodynamic chord, ft 1.2 



TABLE 2 


AIRCRAFT MASS CHARACTERISTICS 


Full Fuel Condition, Without Pilot and Co-pilot 

Weight, lb 3355. 

CG, percent of reference chord 15.9 

Ijj, slug • ft^ 2601. 

ly, slug • ft^ 2052. 

Ij, slug • ft^ 4713, 

Empty Fuel Condition, Without Pilot and Co-pilot 

Weight, lb 2992. 

CG, percent of reference chord 15.7 

I^, slug • ft^ 2526. 

ly, slug • ft^ 2038. 

I^, slug • ft^ 


4608. 


17 


TABLE 3 

INSTRUMENT LOCATIONS RELATIVE TO 
REFERENCE CENTER OF GRAVITY 


Instrument 

Distance forward 
of reference CG, 
ft 

Distance right of 
aircraft centerline, 
ft 

Distance below 
reference CG, 
ft 

a 

3.48 

17,06 

. 833 

B 

3.83 

17,625 

1.4 


-.158 

.885 

.862 


-.158 

.885 

.862 


-.158 

.885 

.862 


158 


.885 


862 
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TABLE 4 

RECORDED PARAMETERS 


Parameter 

Angle of attack, deg 

Angle of sideslip, deg 

Airspeed, ft/ sec 

Altitude, ft 

Pitch rate, deg/ sec 

Roll rate, deg/ sec 

Yaw rate, deg/sec 

Pitch attitude, deg 

Roll attitude, deg 

Yaw attitude, deg 

Normal acceleration, g 

Longitudinal acceleration, g 

Lateral acceleration, g 

Aileron position, total, deg 

Rudder position, deg 

Stabilator position, deg 

Throttle position, left, units 

Engine manifold pressure, right 

Engine RPM, right and left 


Resolution 

.03 

.05 

.15 

.2 

.12 

.12 

.08 

.34 

.39 

.34 

.01 

.002 

.001 

.04 

.06 

.03 

.009 

and left, PSIA .04 


.09 


FULL SCALE WIND TUNNEL DERIVATIVE ESTIMATES 
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FIRST DRYDEN FLIGHT DERIVED DERIVATIVE ESTIMATES 
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APPENDIX A 

MMLE3 AIRCRAFT EQUATIONS OF MOTION 

The following is a list of the equations of motion used in 
the computer program MMLE3, covered in detail in reference 5. 

The longitudinal state, control, observation and extra signal vectors 
are 

X = (a q 0) 
u = (5g) 

extra = (q 6 p r h V) 



The longitudinal observation equations are 
am = Kct (a - ^ q) 

^m = ^ 
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. qS - "an . 

' mg ‘-N — q 

The expansions of the longitudinal force and moment coefficients are 
^ a + Cn,^ ^ + Cm^e *^6 "" ^hiq 

% “ a •" 5e + 

= % 


The approximation of = C,^ is good for low angles of attack. 



APPENDIX B 


DETERMINATION OF MOMENTS OF INERTIA 

The moments of inertia of the PA-30 were experimentally 
determined by the spring oscillation method. The moments of inertia 
(I^, ly, 1 2 ) were determined for empty and full fuel conditions with 
no crew, and landing gear retracted. The vertical center of gravity 
position was also determined from these experiments. 


Vertical CG 

The vertical center of gravity was determined by the single- 
point suspension method (reference 7) . In essence, this method applies 
known loads (Fg^) to displace the suspended aircraft an angle (a) from a 
horizontal attitude. The following equation from reference 7 was used 
to determine the vertical CG position Z from the pivot support: 


F X 
a r a 

tan (a) 


= — C 

W ^ 




Eq. 1 


where is the vertical distance from the supporting pivot axis to the 
aircraft reference line (waterline zero) , The parameters on the right 
side of the equation were measured directly from the test setup. The 
displaced angle, ct, was measured with an inclinometer affixed to the 
aircraft. The test configuration is shown in figure 8. The weight 
included all of the test equipment such as the cradle and supporting 
harness. The test equipment contributions were subtracted from the 
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experimental results. The vertical CG position was found to be about 
4.33 inches below the aircraft's zero waterline reference. 

Moments of Inertia 

The moments of inertia were obtained experimentally by testing 
the aircraft and supporting cradle as a single unit. The equations 
used to solve for the moments of inertia were derived by summing the 
moments about p in the simple model of a spring weight system shown in 
figure 9. The applicable equation listed below was derived in reference 
2 and is a general second order differential equation for damped har- 
monic motion. 


e + 


C . , a^k - Wh ^ 

f e . C — I ) 9 = 0 

^o -^o 


Eq. 2 


where c is a viscous retarding moment and K is the spring constant. 

Equation 2 is a linear differential equation with constant 
coefficients and is solved by the substitution method in reference 2. 
The solution of equation 2 is: 


I 


o 


a^k - Wh 




T 


Eq . 3 


where cOfj is the natural frequency of the system. The term Iq is the 
moment of inertia about the axis of rotation. The natural frequency 
was determined from the time history plots of angular velocity in the 
following manner. The damped frequency, and the damping ratio, C, 
were computed as follows: 


ojjj = 2Trf 
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= 


/4ir^ + <S2 


where f, the frequency of oscillation (cycles/second) , is taken from 
the time history plots, and <S, the logarithmic decrement, is calculated 
using 


1 ^f> 

6 = i £n — 
n Xn 


with n the number of cycles and x the magnitude of the angular velocity 
recorded on the time history plots. 

Then the natural frequency was determined using the equation 

below. 




wd 


Cl - 




The moment of inertia, about the rotation axes, Iq, was then 
transferred to the aircraft CG position using the parallel axis theorem 
to get Ij,g. 


ka^ - Wh w 

ira ~ 5 _ — d^ - 

eg g te 


Eq. 4 


The equipment used in this test included the 256.5 pound 
supporting cradle, the springs, safety cables, and cradle harnesses. 
Anything that moved with the aircraft was taken into account. Only 
half the weight of items affixed to stationary supports was con- 
sidered. Generally most items were small enough to consider their 


so 

moment of inertia about their own center of gravity to be zero. Thus 
the moment of inertia for test equipment, I^e* was simply 

Ite = C'ite) ^ 

where is the weight of the item and is the distance from the 
center of gravity of the item to the oscillation axis. 

The cradle's moment of inertia about its own CG, was large 

enough to be considered, and was found by using the spring oscillation 
method 

Ite (cradle) = (dte)^ 

Equation 4 is used for both the pitch and roll axis tests. The 
pitch test for 1^., pictured in figure 10, used one spring attached to 
the tail tiedown with the center of gravity of the aircraft located 
forward of the pivot point. The roll test for I^, figure 11, used two 
springs affixed, one each, to the wing tiedowns and pivoted directly 
below the center of gravity position. 

The single support cable, or pivot axis, for the yaw inertia 
test configuration, shown in figure 12, was directly above the center 
of gravity position of the aircraft so that the horizontal reference 
line was parallel to the floor and the single supporting cable coin- 
cided with the oscillation axis. This alignment results in zero values 
for distances d and h in equation 4. Thus the equation for yaw is 


h 



'■te 
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The yaw inertia test used four springs horizontally attached with two 
at each wing tiedown point. The yaw analysis assumed the supporting 
cable to be free of torsional moments. 

The spring oscillation method for determining the moments of 
inertia is relatively easy to use. However, many reference points 
exist for measuring distances. So it is a distinct advantage to sim- 
plify the procedure by recording all the vertical and horizontal 
distances to a reference point on the floor and knife edges. 

The tests to determine the moments of inertia were made with 
the landing gear retracted at full and empty fuel conditions. This 
made it possible to interpolate the moments of inertia for different 
center of gravity positions. The moments of inertia are listed in 
Table 2. 
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